International Journal of Advanced Engineering Research and Science (IJAERS) 
https://dx.doi.ora/10.22161/iiaers.6.6.68 


[Vol-S, Issue-6, June- 2019] 
ISSN: 2349-6495(P) 12456-1908(0) 


Study on Wave Type at Water Wave Surface 
Equation Obtained from Kinematic Free Surface 
Boundary Condition (KFSBC) 

Syawaluddin Hutahaean 

Ocean Engineering Program, Faculty of Civil and Environmental Engineering,-Bandung Institute of Technology (ITB), Bandung 40132, 
Indonesia 


Abstract — In this research, water wave surface equation was formulated by integrating kinematic free surface 
boundary condition against time. Then, a study was conducted on the type of wave produced by the water wave 
surface equation. 
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I. INTRODUCTION 

Water wave surface equation of Airy water wave theory 
was formulated by working on small amplitude wave 
assumption (Dean (1991)) that produces sinusoidal wave 
which can be defined that wave height is two times the 
amplitude. Water wave surface formulation was done by 
working on BemouUi equation on water wave surface with 
small amplitude assumption, so the water wave surface 
coincides with stiU water level. Then, KFSBC was done to 
formulate dispersion equation. Therefore, BemouUi 
equation along with KFSBC can produce sinusoidal wave 
type. 

In this research, water wave surface equation was obtained 
by integrating KFSBC against time without working on 
smaU ampUtude wave assumption. In that water wave 
surface equation there are two wave constants that must be 
known their values, i.e. wave constant Gand wave number k 
that are found in the velocity potential solution of Laplace 
equation. The equation to calculate those two variables was 
formulated using KFSBC and surface momentum equation. 
The characteristic of the produced water wave surface was 
studied using Wilson criteria (1963). Based on the criteria 
there are 4 (four) types of wave beginning with wave with 
the smallest amplitude, i.e. Airy’s waves, Stake’s waves, 
Cnoidal waves and Solitary waves. In the deep water, the 
water wave surface equation that was obtained, produces 3 
(three) types of the first wave, depending on the wave 
amplitude as the input. 


II. THE FORMULATION OF WATER WAVE 
SURFACE EQUATION 

Water wave surface equation will be formulated using 
velocity potential equation and KFSBC. Velocity potential 
equation of Dean (1991) is, 

0ix,z,t) = Gcoskxcoshkih + zlsinat .(1) 

xis horizontal axis, zis vertical axis where z = Oat the 
surface of stiU water level, t time, G wave constant, k wave 

2ji 

number, cr = — angular frequency, T wave period and h 


StiU water depth. The equation was obtained by completing 
Laplace equation with variable separation method, where 


coshkih + z)is just a zfunction, so that 
afc(ti+z) 


Rcoshk (ti+z) 

at 


sinhk ih + z) 
ak(h+z) _ Q 


at 


at 

...( 2 ) 


0, thereby 


For aU z value. Forz = r], whererj = riix, t)is the elevation 
of water surface with respect to stiU water level, 
afc ih+T]) 


0 


•(3) 


Equations (2) and (3) are caUed wave number conservation 
equation against time-t. 

With the velocity potential, the equations of particle 
velocity horizontal-xdirection and vertical-z direction can 
be obtained, sequentially as foUows 


u = -= Gksinkxcoshk Qi z)sinat ....(4) 

ax 

w = — — = —Gkcoskxsinhkih z)sinat ...(5) 
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Particle velocity at water wave surface was obtained by 
substituting zwithtj = tjCx, t), whererjis the elevation of 
water wave surface against still water level. 


= Gksinkxcoshkih + r])sincFt ....(6) 

= —GkcoskxsinhkOi + r])sinat .(7) 


Water wave surface equation was formulated by integrating 
KFSBC, y — = — u„— against time-t, . Substitute (6) 

and (7), 

y — = -GksinhkOi + r])coskxsinat 
3t 

— Gkcoshkih + rj)sinkxsinat ^ .(8) 

yis weighting coefficient at weighted total acceleration 
equation with a value around 2.202-3.0. This research uses 
y = 2.483 (Hutahaean (2019a-b)). 


To make the writing easier, the following equations were 
defined 

(rj) = sinhk {h rj) .(9) 

pir]) = coshk Oi -f rj) .(10) 

Substitute (9) and (10) to (8) and by multiplying the right 
side with then (8) becomes, 
a?7 Gk 

— =- p^{ij)acoskxsinat 

— — p (jj) (TsinkxsiTKJt — .( 11 ) 

ycT a^: 

(11) was integrated against time-tby bearing in mind (3), 

i.e. atz = ^ = Oapplies, which means that 

at 

P^iij)andpiij)have constant values against time t,thus the 
integration of the first term right side of the equation was 
completed by integrating sinat, 

Gk 

rj{x,t) = — p^(j]) coskxcosat 
ya 

Gk f dij 

- Phj) asinkx sinat — dt 

ya J 5x 

The integration of the second term right side of the equation 
wiU be completed using partial integration method, as 
follows 


Assume a function off = cosat —. This coscrtfunction was 
' ax 

used so when it was differentiated against time t, sinat^ 
will be formed, i.e. 

at 


5^Tj 


5ij 

asinat -1- cosat - 

5x 5t5x 


This differential equation was multiplied with dtand 
integrated against time t. 


/ 


f ■ . f 

-a J sinat — dt -h J 


cosat - dt 

5t5x 


Substitute /and the first term right side was moved to the 
left and /was moved to the right and both equations were 
divided by cr, 

r 5r] 1 5r] 1 r 

sinat—dt = - cosat -1— cosat - dt 

J 3x a 3x a J 3t3x 

The integration of the second term right side of the equation 

can be completed the same way, but with an assumption 


that 


at^a* 


is a very small number, the integration can be 


completed by integrating just thecosat element. 


/ 


3r] 1 3r] 1 3^r] 

sinat—dt = - cosat -1 —-sinat - 

3x a 3x 0-2 




Substitute the result of integration, 

Gk Gk 

ij{x,t) = — p^{j])coskxcosat - pirj) sinkxsinat 

ya ya 

3ri 1 

— cosat -1— sinat — , 

3x a 3t3xJ 

Working on an assumption that -sinat a very small 
number and can be ignored, 

TjKx.t) = — p^hj) coskxcosat 
ya 

Gk 3ij 

H- p\j])sinkxsinatcosat — 

ycr 3x 

The equation was differentiated against horizontal-x axis. 

3ij Gk 

— =- p, {i])ksinkxcosat 

3x ya 

Gk 3ri 

H- phj) kcoskxsinatcosat — 

ycr 3x 

Gk ^ ^ 3^ri 

H- Pin) sinkxsinatcosat —- 

-r- / 


ya 


3x^ 


or 


^ — Pi(T])ksinkx ■+ — p(j])kcoskxsinat^ 

3 ^ 77 \ 

P {rj)sinkxsinat cosat .(12) 


ax 

— 

ya 

In accordance with the provision at the velocity potential 
equation where there is tfunction only,xfunction 
onlyandzfunctiononly, then at the water wave surface using 
variable from velocity potential equation, water surface 
equation also has variable that is tfunction only anxfunction 
only, so it can be stated that the general form of (12) is, 

3t] 

3x 

In this equation /(x)is just a function of x. /(x)in (12) is. 


/ (x) cosat 
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/ Gk G k 

fix) = I- fJ^iii)ksinkx H- pirf) kcoskxcosatf ix) 

\ ya ya 

ya 3%/ 

Bearing in mind (3) fJ^ir]) and fir]) should be constant 
numbers against time tand against horizontal x axis. In the 
second and third terms right side of the equation there is a 
function of time t, i.e.cosat , then the term shouldn’t be 
there. Therefore/ (x) is, 

Gk 

fix) = - f^iri)ksinkx 

ya 

Thereby 
drj Gk 

— =- p, irj) ksinkxcosat 

3x ya 

Substitute this equation to (11), 

5r] Gk 

— =- p^lj]) acoskxsinat 


+ 


(—) pirj)p^ir])aksin^kxsinatcosat 
\ya > 


is integrated against time t. Integration will be completed 

by bearing in mind wave number conservation equation (3), 

as in the previous section, thus integration was done by 

integrating only sinusoidal element and bearing in mind that 

sinatcosat =-sin2at , 

2 

Gk 

r]ix,t) = — pAr]) coskxcosat 
ya 

1 fGkA 

——j pirpPArpksiiAkxcosZat 

..(13) 

To determine the values of /iC??), pir]), an approach was 
done that water surface equation is sinusoidal, i.e. 
rjgixit) = Acoskxcosat .(14) 

Therefore, the approach of hyperbolic function is Piirf) = 

l^ii^oX l^i^') = PiTJo)- 

riix,t) = —PM coskxcosat 
ya 

IfGkA 

——j p{r]Q)P^{r}Q)ksin^kxcos2at 

-(15) 


ya ^ 3x 

At the characteristic point, i.e. at coskx = sinkx = 
cosat = sinat, the equation becomes, 

^ = .<>«> 

Surface momentum equation where convective acceleration 
term is ignored is. 


au„ 


-g- 




..(17) 


At the characteristic point, where "^was obtained from (6) 
equation yaGp = ^riwas obtained 


or, 

G = 


gA 


..(18) 


y<yp{^) 

Substitute (17) to (15), 

y^a^ = gktanhk (h -l-^ — 


..(19) 


was obtained. In the deep water, tanhk (h + jj = 1, then 
dispersion equation in the deep water is, 

gk^A 


y^a^ = gk-'- 


( 20 ) 


This equation is quadratic equation of wave number 
fcwherefccan be calculated using simple method, i.e. finding 
the square root. Therefore, the calculation of G and fcwas 
done by calculating fcwith (20) and then Gwas calculated 
with (18). 

As an illustration of the result of the calculation, wave 
lenght atLwas calculated whereL = ^ , for wave with wave 

period of 8 sec., wave amphtude 0.6 m, 0.8 m and 1.0 m 
The result of calculation, shown on Fig. Ishows that there is 
an impact of wave amplitude on wavelength, i.e. the bigger 
the wave amplitude, the shorter the wavelength. 
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Water depth h (m) 


III. THE CALCULATION OF GANDfc 

To calculate Gandfc, two equations are needed, i.e. KFSBC 
and surface momentum equation. 

By working on (14) on (11), 


at 


- — pM(^coskxsinat 
ya 


A=0.60 m A=0.80 m 
-A=1.00 m 

Fig.l. The impact of wave amplitude A on wavelength L 

2k 
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In (20), there is a maximum value of wave amplitude, i.e. 
when determinant value D = “ 4(y^cr^) = 0. In 

this condition wave amplitude A = , with the value, 

.( 21 ) 


At wave amplitude maximum where determinant = 0 , 
wave number becomes 
feo = 7^ .(22) 

^max 

In the deep watertan/tfcg (^hg +^) = tanhkghg (^1 + 
—) » tanhkghg = 1. Consider— « 1. 
wheretanh (I.OBtt) = 0.999937. Then deep water depth hg 
ho =^= I.OSttTI^,, .(23) 

The use of coefficient 1.65 in (23) is the result of breaker 
depth calibration in shoaUng-breaking model against 
breaker-depth from SPM(1984). The model is not discussed 
here. 



Fig.2. Graph of the value ofGas a function of wave 
amplitude A 


result shown on Fig.l. should be read as a wave with wave 
constant values Gand determined wave amplitude. 
Therefore, Fig.l. is not the value of wave amplitude 
kmoving from deep water to shallow water. 

IV. TYPE OF WAVE PRODUCED BY WATER 
WAVE SURFACE EQUATION 
4.1. General Shape 

The general shape of water wave surface produced by(15) is 
shown on fig.3. The shape of water wave surface is not 
symmetrical between wave crest and wave trough where 
hmax > IhminI’ Vmax^^ '^^e elevation of wave crest above 
stiU water level (line of elevation z = 0) , whereas is 
the elevation of wave trough below line z = 0. Wave 
profile like this is called cnoidal wave profile 
(Wilson(1963)). With wave height H = g^ax ~ hmin ■ 



Fig.3. General shape of water wave surface. 


Fig.2 shows graph of value Gas the function of wave 
amplitude, for wave with wave period of 8 sec in the deep 
water. It shows that the value of Ggrows bigger as the wave 
amplitude grows bigger, but there is a wave amplitude value 
where Greaches maximum value, i.e. the value of 
Gdecreases as wave amphtude grows bigger, until it has 
negative value at wave amplitude equals toA^^^. There is a 
similar phenomenon for other wave period, whereG^^j-was 
achieved at wave amplitude of 0.91 Therefore, for a 

wave period, the wave amplitude value that should be used 
is 


In (22) and (23) A^^^was used in (24). Bearing in mind (3), 
the calculation of wave number kfor a wave moving from a 
certain depth to another depth cannot be calculated with 
(19), it should be with wave number conservation law. The 


4.2. Wilson Criteria (1963) 

Wilson (1963) classified wave type based on the value of 
. Airy waves or sinusoidal waves types have a value 
of < 0.505. This shows that the shape of wave is 

symmetrical, wheretj^^j. « Itlminl- TyP® Stoke’s wave type 
is StiU quite close with the sinusoidal character. 
Furthermore, cnoidal waves type is not symmetrical at aU 
where there is a form of wave with the wave trough that 
almost coincides with the line z = Oi.e., where«= 1.0. 

H 

This type of wave is called soUtary waves type. 


Table.1: Wilson criteria (1963) 


Wave Type 

Vmax 

H 

Airy waves 

< 0.505 

Stoke’s waves 

< 635 
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Cnoidal waves 

^ Vmax 

0.635 <-< 1 

H 

Solitary waves 

= 1 


Fig.4. is the shape of wave surface sinusoidal for wave 
period 8 sec., which shows that wave crest and wave trough 
are symmetrical, = ItJminl- 


4.3. Sinusoidal Type (Airy Waves Type) 

Table (2) shows the result of the calculation with a very 
small wave amplitude where the value of = 0.503was 
achieved, which means that the wave profile is not really 
symmetrical, wave crest is stiU bigger than wave trough, but 
it is very close with the symmetry, where the value of “ = 
2was achieved, oxH = 2A . Therefore, it can be determined 
that as a criteria Airy wave type or sinusoidal wave type is 
H = 2A. The example of sinusoidal wave profile is on 
Fig.4., wherewave crest with wave trough is symmetrical, 
i.e. rjjnax ~ ^Vmin I- 


Table.2: Wave amplitude for sinusoidal wave, at deep 
water. 


T 

A 

Vmax 

H 

H 

(sec.) 

(m) 

H 

(m) 

A 

8 

0,017 

0,503 

0,035 

2 

9 

0,02 

0,503 

0,041 

2 

10 

0,023 

0,503 

0,047 

2 

11 

0,027 

0,503 

0,053 

2 

12 

0,03 

0,503 

0,06 

2 

13 

0,033 

0,502 

0,067 

2 

14 

0,037 

0,502 

0,073 

2 

15 

0,04 

0,502 

0,081 

2 


0.02 



Fig.4. Sinusoidal wave T = 8 sec., A = 0.017 m, at deep 
water. 


4.4. Cnoidal Type 

Table (3) and Table (4) show the result of the calculation of 
wave profile characteristic with wave amplitude 0.523^4^^^ 
and . 

With wave amplitude of0.523y4^a;c, ’^'"°^ = 0.635was 

H 

obtained, which, according to Wilson criteria is the 
maximum limit of Stoke’s wave type. The value of - = 
2.099, is stUl quite close with 2. 

Table.3: Wave profile at wave amplitude A = 0.523^4^^^^ 


T 

A 

Vmax 

H 

H 

(sec.) 

(m) 

H 

(m) 

A 

^max 

7 

0,47 

0,635 

0,986 

2,099 

8 

0,614 

0,635 

1,288 

2,099 

9 

0,777 

0,635 

1,631 

2,099 

10 

0,959 

0,635 

2,013 

2,099 

11 

1,161 

0,635 

2,436 

2,099 

12 

1,381 

0,635 

2,899 

2,099 

13 

1,621 

0,635 

3,402 

2,099 

14 

1,88 

0,635 

3,946 

2,099 

15 

2,158 

0,635 

4,529 

2,099 


With wave amplitude ofA^^^, = 0.803 was obtained, 

H 

with the value of^ = 2.726 > 2with a deviation of 36.3 % 

A 

of 2. 


Table.4: Wave profile at wave amplitude A = A^^^ 


T 

A 

Vmax 

H 

H 

(sec.) 

(m) 

H 

(m) 

A 

7 

0,899 

0,803 

2,449 

2,726 

8 

1,174 

0,803 

3,199 

2,726 

9 

1,485 

0,803 

4,049 

2,726 

10 

1,834 

0,803 

4,998 

2,726 

11 

2,219 

0,803 

6,048 

2,726 

12 

2,641 

0,803 

7,198 

2,726 

13 

3,099 

0,803 

8,447 

2,726 

14 

3,594 

0,803 

9,797 

2,726 

15 

4,126 

0,803 

11,246 

2,726 


The profile of wave with wave period 8 sec., wave 
amplitude A = 1.174 m or^l can be seen on 
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Fig.5. The cnoidal profile that was formed is very clear and 
almost perfect. 



Fig.5. Cnoidal wave T = 8 sec., A = 1.174 m, at deep 
water. 


V. CONCLUSION 

In the deep water, there are three types of wave produced by 
kinematic free surface boundary condition, i.e.Airy waves, 
Stoke’s waves dan cnoidal waves. At a very small wave 
amplitude. Airy waves type was produced with symmetrical 
wave profile between wave crest and wave trough and wave 
height has a value twice of wave amplitude value. By 
enlarging wave amplitude, Stoke’s wave type will be 
obtained, where the wave height is stiU close to twice the 
wave amplitude. Cnoidal wavetype was obtained by 
enlarging wave amplitude more than the wave amplitude at 
Stoke’s waves, where at cnoidal waves the value of wave 
height is bigger than twice the wave amplitude. 

The second conclusion is a calculation with an assumption 
that wave height is twice wave amplitude for a large wave 
amplitude is inaccurate 

As has been stated, the research was done at deep water. 
Therefore, the next research should be on wave profile at 
shallow water, particularly at breaker point. This is 
important to interpret the result of shoaling breaker analysis 
with wave amplitude variable. It is also to examine wave 
height produced by a model using wave height as its 
variable by defining wave height that is twice wave 
amplitude 
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